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and the Nature of the Self- Trapping Transition 
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We analyze electron-phonon correlation functions mea- 
sured in ID polaron ground states of the Holstein Hamiltonian 
using the Global-Local variational method. The spatial col- 
lapse of electron-phonon correlations is found to occur in con- 
cert with transition behavior in other polaron properties, pro- 
viding mutually confirming evidence for a self-trapping line in 
ID. The spatial extent of electron-phonon correlations is used 
to quantify polaron size, and is analyzed over a wide range of 
parameters. Distinct scaling behaviors are found to be char- 
acteristic of the region below the self-trapping transition and 
above it, contrary to some widely-held expectations and lead- 
ing naturally to the notion of the polaron size as an order 
parameter for a self-trapping transition that becomes critical 
in the adiabatic limit. 

PACS numbers: 71.38.+i, 71.15.-m, 71.35.Aa, 72.90.+y 



Polarons are the quantum quasiparticles describing the 
states of excitations correlated with the deformation or 
polarization quanta of a host medium. Perhaps the sin- 
gle most intuitively-accessible polaron property is its size; 
however, the size of the polaron is one of its properties 
that is most elusive. Generally, we do not measure it di- 
rectly, but typically impute to it some qualitative char- 
acter (e.g., "large" or "small") by interpreting measure- 
ments of more directly-accessible properties within some 
theoretical view of the problem, often reduced to rules of 
thumb. 

In this paper, rather than rely upon rules of thumb 
and qualitative characterization, we subject the notion 
of polaron size to direct quantitative test. A large vol- 
ume of numerical data is analyzed to put numerical value 
to measures of the polaron width over a broad range of 
the polaron phase diagram, and the essential characteris- 
tics of polaron size are mapped out. The volume, scope, 
and accuracy of our data makes it possible to discern 
trends that indicate new "rules of thumb" that challenge 
existing thinking about polaron size. 
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The size of the polaron is not uniquely-definable owing 
both to the inherent complexity of the polaron quasipar- 
ticle and to the strictly delocalized character of polaron 
eigenstates in translationally-invariant lattices. There is 
nonetheless no shortage of opinion regarding the spatial 
structure of the polaron, much of which is tied to the 
adiabatic approximation ; some of this is based on 
localized-states, the balance on delocalized states built 
up from localized states. 

We work in a strictly delocalized-state picture, wherein 
we construct polaron Bloch states as superpositions of 
localized functions ( "form factors" ) that are the vehicles 
through which local electron-phonon correlation struc- 
ture is conveyed into the delocalized momentum eigen- 
state. The quality of result is directly attributable to 
the appropriateness of the local form factor. Much work 
has been done on polaron states as localized states, from 
which point of view the form factors we use are them- 
selves viable localized polaron states. While some simi- 
larity of result between such approaches can be expected, 
it is essential to the following that such similarity is lim- 
ited. When integrated into the delocalized Bloch state, 
the "identity" of the form factor is lost and one can no 
longer characterize the properties of the quantum state 
in terms uniquely associated with it. Local properties 
can be probed, however, with correlation functions that 
reveal the internal structure implicit in the delocalized 
state; thus, we characterize polaron size through the use 
of electron-phonon correlation functions that reveal the 
structure of the "phonon cloud" or co-moving lattice dis- 
tortion that surrounds the electron at every instant. 

We use the Global-Local variational method to deter- 
mine a large number of such ground states over broad 
regions of the polaron phase diagram. The particulars 
of the Global-Local method and other results obtained 
through its use have been detailed elsewhere ]7|-|l4[] . The 
Global-Local method has been demonstrated to be con- 
sistent with both weak-coupling and stong-coupling per- 
turbation theories in the appropriate regimes, and over 
the regimes studied here, to be in broad quantitative 
agreement with other high-quality methods such as clus- 
ter diagonalization [fl5|-ll8| , density matrix renormaliza- 
tion group |l9| , and Monte Carlo simulations [^,|T) . 

We limit ourselves to ID and use the traditional Hol- 
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stein Hamiltonian M denned by 

H = - J y^g^(g n+1 + a n _i) + huj blb n 

n n 

(bl + K) , (1) 

n 

in which a' n creates a single electronic excitation in the 
rigid-lattice Wannier state at site n, and 6^ creates a 
quantum of vibrational energy Hlo in the Einstein os- 
cillator at site n. The hopping matrix element con- 
necting nearest neighbors is given by J, and g is the 
electron-phonon coupling strength. We will use Tiuj as 
the unit of energy, in terms of which the two independent 
model parameters are defined (hopping integral J/hui 
and electron-phonon coupling g). 

The substance of our results should not depend greatly 
on precisely what correlation function we use, provided 
that electron-phonon relationships are spatially resolved. 
We use the function 

Cr = Y Q E(°J.°n( 6 n+r + ' ^ 

" n 

normalized such that ^ r C r — 1. Examples of C r at gen- 
eral J/hiv and g obtained by the Global-Local variational 
method are shown in Figure HQ. These correlation 
functions are strongly exponential over four decades of 
amplitude in the site index r except in the intermediate 
coupling region where mild deviations from exponential 
decay are found. 



The dependence of individual C r on model parameters 
can be used to determine the location of the self-trapping 
transition. Through a curvature analysis of a large vol- 
ume of Global-Local data, we have identified the self- 
trapping transition with the point of most rapid decrease 
in C' r w.r.t. g at fixed J/Tilo for several sites near the 
center of the polaron. Figure |^ shows that these loci for 
CV=o,i.2,3 cluster tightly about the same self-trapping line 
(ffsr) as has been previously identified through numer- 
ous other polaron properties [p|, ^0|Jl^ , [l^ ] , assuring that 
all these analyses are identifying the same fundamental 
event. The apparent self-trapping trends in all these are 
well located by the empirical curve gsT = 1 + \/ J/huj as 
shown in Figure |^. 
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FIG. 1. Electron-phonon correlation function CV as deter- 
mined by the Global-Local method. J/Huj = 7.0 for assorted 
values of g. Open symbols lie in the small polaron region, 
solid symbols in the large polaron region. 



FIG. 2. Polaron phase diagram, showing the location 
of the self trapping line gsr and the specific loci due to 
Global-Local measurements of C r =o, 1,2,3- 

The self-trapping line gsr is a very important feature 
in all these considerations because it sets the coupling 
scale around which the most physically-significant fea- 
tures are organized. At large g, above the self-trapping 
transition, the polaron is extremely compact up to the 
self-trapping line where a fairly rapid spreading of corre- 
lations commences; however, the spreading of the corre- 
lations saturates at a finite width in the weak-coupling 
limit. This saturation phenomenon proves to be a crucial 
finding at variance with many widely-held views about 
the nature of large polarons j8| Jll|jr^ , |l^] . 

Our basic tool for quantifying polaron size is the vari- 
ance of the correlation function 



<t 2 = 5> 2 CV 



(3) 



An important first result from weak-coupling pertur- 
bation theory @,p,0,ra is 
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t 2 2 

lim a = <r n = - — 

5^0 ftW 



(4) 



This limiting behavior is what is seen in Figure |], where 
the spreading of electron-phonon correlations saturates 
at a finite width as coupling vanishes. 

On the other hand, above the self-trapping transition 
it can be shown (for example, using strong-coupling per- 
turbation theory HQ) that 



lim 

J/huj — >oo CJq 



= 



9 > 9ST 



(5) 



This result involves the fact that being "above" the self- 
trapping transition while J/hw diverges requires that g 
diverge as well in a sufficiently rapid manner. 

With these properties in mind, we scale our variance 
data so that they can be meaningfully organized relative 
to the self-trapping transtion. The weak-coupling vari- 
ance <7q appears to be the natural quantity with which to 
scale variances, and gsr the natural quantity with which 
to scale the electron-phonon coupling. 

Using the Global-Local variational method, we have 
computed a 2 for approximately 1200 polaron ground 
states and have compiled these results in Figure [| us- 
ing the above scaling conventions. 
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FIG. 3. Scaled correlation variance a 2 l<j\ as results from 
the Global-Local variation. Variance data are truncated as 
discussed in the text. 

These data deteriorate at the weak coupling end of 
each curve where the Global-Local method encounters 
systematic difficulties. The variance as a measure of spa- 
tial structure amplifies contributions from tails of the po- 
laron wave function, indiscriminately amplifying any er- 
rors lurking in that part of the numerical solution. This 
poses problems for variational methods, since the varia- 
tional energy is least sensitive to the detail of such tails, 



with the result that relative errors are largest in this part 
of the solutions. Similarly, any finite-size effects present 
are most deleterious in the wave function tails. For such 
reasons, only the descending part of each numerically- 
determined curve is to be regarded as reliable; we have 
left a small ascending arc at the weak-coupling end of 
each curve only to render clearly where and how the de- 
terioration of results sets in. 

Focussing therefore on the descending data in Figure |^, 
we are led to distinguish two apparent facts: I) that the 
variances above the self-trapping transition (g > gsr) 
appear to be converging toward zero (consistent with 
SCPT), and 2) that the variances below the self-trapping 
transition (g < gsr) appear to be converging toward fi- 
nite values (consistent with WCPT). These findings in- 
dicate a difference in scaling properties above and be- 
low the transition that suggests the emergence of critical 
behavior in the adiabatic limit. Relative to such emer- 
gent critical behavior, the limiting variance curve below 
the transition has the characteristics of an order param- 
eter. What is most significant in this is not the detailed 
form of that limiting curve, which we do not suggest our 
data is sufficiently complete or accurate to determine, 
but merely that such a limiting curve appears to exist. 

A fair characterization of the apparent trends in con- 
vergence is contained in the function 



\9stJ 



= 4/1- 



9st 



9 < gsr 



9 > gsr 



(6) 



as indicated by the beaded guide curve in Figure || and 
under a different scaling convention (see below) in Fig- 




FIG. 4. Order parameters characterizing the adiabatic 
limit. Solid curve, Eqn. [| Dashed curve, Eqn. ^. 
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Thus, we find that the polaron size below the self- 
trapping transition is roughly proportional to that found 
in the weak coupling limit, with the constant of propor- 
tionality dependent primarly on the "distance" from the 
self-trapping transition. 

This finding is consistent with and complementary to a 
well-known result of adiabatic strong-coupling theory |EJ 
that identifies a soft vibrational mode whose existence is 
synonymous with the existence of the self-trapped state. 
The soft mode frequency uj is interpretable as an order 
parameter for the self-trapping transition. In terms of the 
composite parameter A = g 2 huj/2J common in strong- 
coupling theory, 
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= \ 1- 



x 2 
A 2 



A < Xst 



A > Xst 



(7) 



with X ST = 1/2 in ID. (See Figure |.) 

It is perhaps not surprising that the polaron size should 
have the characteristics of an order parameter, since 
the essential physical notion underlying the self-trapping 
transition is that of spatial collapse; however, accord- 
ing to some widely-held and long-established views, such 
collapse is not expected to occur in ID. Under the adia- 
batic approximation in ID, the size-scaling behavior that 
we have here found to be limited to the region above the 
transition is found instead at all coupling strengths, with- 
out undergoing any transmutation into a distinct weak- 
coupling behavior. Consequently, the "rules of thumb" 
given by the adiabatic approximation to characterize the 
size of the large polaron in ID disagree with the "rule of 
thumb" that can be taken from Eqn. ||. 

The discrepancies between the traditional notion of 
the large polaron and what we find here are perva- 
sive and, in regimes, severe. Necessarily, we are led to 
supplant the traditional notion of the large polaron as 
a broad "soft" state that responds rather strongly to 
changes in electron-phonon coupling strength with that 
of a less broad, "stiff" state that is relatively insensitive 
to changes in electron-phonon coupling strength. 

Although the picture of the self-trapping transition 
that emerges from this work has so far been demonstrated 
only in ID, many of the elements and relationships con- 
tributing to it arc already known to hold essentially un- 
changed in higher dimensions Jn],[u|^0|,^l],|3) . We expect 
<7 2 /<7q to continue to constitute an order parameter in the 
adiabatic limit and for the order parameter to provide a 
reasonable "rule of thumb" for the correct size scaling of 
the large polaron in each dimension. 
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